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ABSTRACT 

We present measurements of reddening due to dust using the colors of stars in the Sloan Digital Sky 
Survey (SDSS). We measure the color of main sequence turn-off stars by finding the "blue tip" of the 
stellar locus: the prominent blue edge in the distribution of stellar colors. The method is sensitive to 
color changes of order 18, 12, 7, and 8 mmag of reddening in the colors u — g, g — r, r — z, and i — z, 
respectively, in regions measuring 90' by 14'. We present maps of the blue tip colors in each of these 
bands over the entire SDSS footprint, includin g the new d u sty so uthe rn Galactic cap data p rovided by 
the SDSS-III. The results disfav or the bes t fit (T Donnelll (Jl99l and lCardelli"eTa/] (Ella) reddening 
l aws, but are well described by a iFitzpatrickl (jl999h reddening law with Ry = 3.1. The lSchlegel et all 
(1998|, SFD) dust map is found to trace the dust well, but overestimates reddening by factors of 1.4, 
1.0, 1.2, and 1.4 in u — g, g — r, r — z, and i — z, largely due to the adopted reddening law. In select 
dusty regions of the sky, we find evidence for problems in the SFD temperature correction. A dust 
map normalization difference of 15% between the Galactic north and south sky may be due to these 
dust temperature errors. 

Subject headings: dust, extinction — ISM: clouds — Galaxy: stellar content 



1. INTRODUCTION 

Most astronomical observations are affected by Galac- 
tic interstellar dust, whether as a source of foreground 
light in the microwave, far-infrared and gamma-ray wave- 
length regions or as a cause of ext inction in the infrared 
through ultraviolet ([Drain e 2003). Characterizing the 
properties of the dust and accounting for its effects on 
observations is then a central problem in astronomy. 

Dust is formed as stars burn nuclear fuel to heavy ele- 
ments and emit these elements in stellar winds or in more 
violent eruptions, and t hese elements are reprocessed in 
the interstellar medium (jDrain e 2009). The distribution 
of dust is correspondingly correlated with the hydrogen 
in the inters tellar medium (ISM) , though the dust-to-gas 
ratio varies. iBurstein and Heilesl {1978, BH) took advan- 
tage of the correlation between H i and dust column den- 
sity to make the first widely used map of dust column 
density, combining H I emission and galax y counts. 

The BH dust map was superseded by the Schl egel et all 
(1998, SFD) dust map, which took advantage of the 
full-sky far-infrared (FIR) data provided by IRAS and 
DIRBE, which are dominated by thermal emission from 
the Galactic dust at wavelengths of 100 jam and longer. 
After correction for dust temperature using FIR color ra- 
tios, these maps trace dust column density. Calibration 
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of dust column density to color excess E(B — V) was per- 
formed using the colors of a sa mple of 389 galaxies with 
Mg2 indices (jFaber et a/.lll989f ). which were also used to 
test the performance of the resulting map. 

The advent of large astronomical surveys permits 
stronger tests of the SFD dust map and the assumptions 
used to construct it. We test the SFD dust map using 
the color s of stars from the Sloan Digital Sky Survey 
(SDSS) (jYork et al\ 12000|) including the recent ly com- 
pleted SDSS-III imaging ([Weinberg et a/.ll2007D . which 
covers about 2000 deg 2 in the southern Galactic sky. 
This allows us to more tightly constrain the SFD normal- 
ization and the dust extinction spectrum, or "reddening 
law," over the SDSS bands, as well as to provide a map of 
color residuals that point to problems with the SFD dust 
map, and particularly with the temperature correction. 

Previous tests of SFD have usually found that 
SFD overpredicts extinction in high-extinction r e gions . 
Shortly after its introduction JArce and Goodman! (|1999) 
found that SFD extinction was too high by 40% in Tau- 
rus, using star counts, colors, and FIR emission. Like- 
wise, studies using globular cluster photometry, galaxy 
counts, and NIR colors have found that SFD overpre- 
dicts extinction by a similar fraction in other dusty re- 
gions ([Sta nek 1998; Chen et al 1999; Yasuda et a/1l2007l : 
iRowles and Froebrichl 120091 ). iCambresy et all (|2001f ) ex- 
plore the link between the SFD overestimation and dust 
temperatu re using star c o unts in the Polaris Flare. At 
\b\ < 40°, Dobas hi et al\ (|2QQ5[ ) use optical star counts 
to conclude that SFD overpredicts extinction by a fac- 
tor of two or more. We perform similar tests to these 
in regions of lower extinction than had been previously 
possible, taking advantage of the high quality and depth 
of the SDSS st ellar photometry, a nd co mplementing the 
recent work of iPeek and Grave 1 (|2010h . who use SDSS 
galaxy spectra. In these regions, we do not find that 
SFD overpredicts reddening by a large factor; rather, we 
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find that SFD overpredicts reddening by about 14% in 
B — V, though, because of the reddening law adopted by 
SFD, this varies from color to color. 

The colors of stars vary substantially with stellar type, 
and with location in the Galaxy due to the effect of 
metallicity on color. Nevertheless, we find that the colors 
of the most blue main sequence stars in old populations — 
the main sequence turn-off (MSTO) stars — are remark- 
ably stable over the sky, and that we can empirically 
model their slow spatial variation. We therefore present 
measurements of the colors of the "blue tip" of the stel- 
lar locus, and use them to c onstrain the SFD dust map. 
This work is akin to that of lHigh et ali (j2009[ h in which 
the color of stellar populations is also used as a sort of 
color standard. 

In $2] we describe the data sets used in this work: the 
SDSS imaging data and the SFD dust map. In we 
present our method for measuring the blue tip of the stel- 
lar locus, and the corresponding maps of blue tip colors 
on the sky. In £|4j we present fits of the SFD dust map to 
the blue tip colors. In $5] and $6j we discuss these results 
and conclude. The blue tip maps and measurements can 
be found at the web site skymaps . info/bluet ip. 

2. DATA 

2.1. The SDSS 

The SDSS is a digital spectroscopic and photomet- 
ric survey, that, with the additional south Galac- 
tic cap (SGC) data provided by the SDSS-III, covers 
just o yer one third of the sky, mostly at high lati- 
tudes (Abazaiian ~eTal\\200^ ). The SDSS provides near- 
simultan eous imaging in five optical fil ters: u, g, r, z, 
and z (jGunn et a/.lll998HFukugita et a/.lll996D . The pho- 
tometric pipeline has uniformly reduced data for about 
10 8 stars. The SDSS is 95% complete up to magnitudes 
22.1, 22.4, 22.1, 21.2, and 20.3 in u, g, r, i, and z. The 
SDSS imaging is performed in a drift-scanning mode in 
which the shutter is left open while the telescope moves 
at constant angular speed across the sky. The result- 
ing imaging "runs" tend to cover strips of sky that are 
long in right ascension and narrow in declination. These 
runs are divided into fields, which are 13. 5' by 9' in size. 
We use SDSS data that have been photometrically cal- 
ibrated according to the "u bercalibration" procedure of 
iPadmanabhan et all (|2QQ8h . 

2.2. The SFD Dust Map 

The SFD dust map is a map of thermal emission from 
dust based on the IRAS 100 /im maps. The IRAS 100 /im 
map underwent three major processing steps before being 
used in the SFD dust map: it was destriped, zodiacal- 
light subtracted, and calibrated to match COBE/DIRBE 
at degree scales. The first and third of these steps were 
necessary because the zero point of the IRAS 100 /am 
detector varied over time, imprinting the IRAS scan pat- 
tern on the data in the form of stripes, and making the 
overall zero point of the data uncertain. The zodiacal- 
light subtraction is needed to remove the signature of the 
hot, bright, local interplanetary dust from the IRAS and 
DIRBE data, which, though bright, contributes negligi- 
bly to the reddening. The DIRBE 100 /im to 240 /im 
flux ratio is used to constrain the dust temperature and 
to transform the 100 jam flux map to a map propor- 
tional to dust column density. The destriped IRAS data 



have resolution of 6 arcminutes, while the temperature 
correction is smoothed to degree scales. The destriped, 
temperature-corrected IRAS 100 jam map is normalized 
to a map of E(B — V) by using a sample of 389 galax- 
ies with M g2 indices and B — V colors, for use as color 
standards (jFaber et ali 119891 ). 

3. THE BLUE TIP OF THE STELLAR LOCUS 

Main sequence stars are confined to a tight locus in 
color-color space, as is apparent from a typical SDSS 
color-color diagram (Fig. [I]). The locus has a sharp blue 
edge at the color of the MSTO stars, blueward of which 
there are only rare blue stragglers, white dwarfs, and 
quasars. The intrinsic color of the MSTO is set by the 
properties of the stellar population: primarily, metallic- 
ity and age. The observed colors of the MSTO are the 
intrinsic colors, altered by reddening due to dust and 
by systematic color shifts due to imperfect calibration of 
the observations. Empirically, the intrinsic color of the 
MSTO is slowly varying in space; the observed color of 
MSTO stars is therefore a good probe of the reddening 
due to dust. 

Population effects cause the color of main sequence 
stars to vary as a function of the magnitude of stars 
probed. At faint magnitudes and high latitudes, metal- 
poor halo stars dominate, rendering the MSTO stars 
bluer than at bright magnitudes and lower latitudes, 
where redder disk stars dominate (Fig. [2]). Moreover, the 
shape of the blue edge changes with magnitude: owing to 
photometric errors, the edge is sharp for bright stars and 
becomes blurred at fainter magnitudes. Accordingly, the 
typical color of MSTO stars depends on magnitude. 

Stars of similar metallicities may clump together in the 
halo, changing the observed color of the blue tip. Because 
we do not have a good way to model such effects, we 
cannot distinguish such clumps from reddening due to 
dust in the blue tip maps. 

3.1. Measuring the Blue Tip 

We perform fits to the location of the blue tip of the 
stellar locus in each SDSS field (or group of SDSS fields). 
To avoid biasing the fits, we select stars in each field for 
analysis in a reddening-independent way, and then find 
the location of the blue tip in the resulting sample. 

Care must be taken when imposing cuts on the stars 
in each field to avoid biasing the resulting blue tip color. 
Because extinction changes the observed magnitudes of 
stars, a flat cut on magnitude would bias the measure- 
ment, as the distance range probed would change with 
the dust column. Instead, we perform cuts on D, the g 
magnitude of a star projected along the reddening vector 
to zero color: 

D = S~T^(9-r). (1) 

Provided A g /(A g — A r ) is accurate, the set of stars in a 
particular range of D is independent of the dust column 
towards the stars. D is related to the distance modu- 
lus of the MSTO star s (though not redder stars). Using 
the Uuric et ali (|2008f ) "bright" photometric parallax re- 
lation, MSTO stars with r - i = 0.1 and g - r = 0.2 
have M g = 5.3. Accounting for the projection from 
g — r = 0.2 to zero color along the reddening vector, 
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Figure 1. Color-color diagrams for point sources from the SDSS, within two degrees of the north Galactic pole. The sharp cutoff in 
stellar density blueward of the "blue tip" is prominent. In u — g, the points blueward of the blue tip are mostly quasars and white dwarfs. 
The blue tip is at about 0.8, 0.2, 0.1, and 0.0 mags in u — g, g — r, r — i, and i — z, respectively. 




Figure 2. Color-magnitude diagrams from the SDSS, for the same stars as in Fig. ^ The location of the blue edge of the stellar locus 
depends on magnitude, especially in the bluer bands, where bright, nearby, more metal-rich stars are redder than more distant halo stars. 
In u — g, blueward of the blue edge lie mostly quasars and white dwarfs. The location of the blue tip is as in Fig. [T] 



D w g — M g + 4.6. In this work, we frequently use the 
range 10 < D < 19. Including the effect of saturation 
on the bright end, this restriction on D selects stars with 
distances between about 1 and 8 kpc. 

With a reddening-independent population of stars in 
hand, we measure the color of the blue edge of these 
stars in each SDSS color (Fig. [3]). The stellar colors are 
modeled as drawn from a probability distribution P(x), 
taken to be a step function convolved with a Gaussian 
to reflect the photometric uncertainty and the intrinsic 
width of the blue edge. The maximum-likelihood loca- 
tion for the step is denoted the "blue tip" of the locus. 
Specifically, the form of the probability distribution is 
taken to be 



1 + erf 



X — Xq 



V2i 



(7 p 



(2) 



where erf is the Gaussian error function, op corresponds 
to the width of the edge and F sets a floor to the prob- 
ability distribution, to render the fit insensitive to the 
occasional white dwarf or blue straggler. In this work we 
use F = 0.05. The variable xo is the only free parameter 
in the fits, and gives the color of the blue tip. The final 
results of this work are insensitive to substantial changes 
in the floor F, from 0.01 to 0.05. The formal statisti- 
cal uncertainty in the blue tip color is computed for each 
measurement by fitting the likelihood near the maximum 
likelihood to a Gaussian. 

The width of the blue edge is set by a combination 
of the intrinsic width of the edge and the typical mea- 
surement error in star colors in the SDSS. We adopt the 
values 0.05, 0.05, 0.025 and 0.025 for the width parame- 
ter crp in the colors u— g, g — r, r — i and i— z, respectively. 



Because the photometric uncertainty becomes larger at 
faint magnitudes, at high extinctions the observed blue 
edge is expected to broaden. Despite this, ap is kept 
constant in the fits, introducing the possibility of bias 
in the measurements. Nevertheless, as we look only at 
relatively bright ranges of D and low extinction, this is a 
minor effect. Simulations with a mock star catalog ( §4.4|) 
verify that for the range 10 < D < 19 used in this work, 
the introduced bias in the recovered ratio of blue tip color 
to E(B — V) is less than 2% in u — g and smaller in the 
other colors. 

The statistical uncertainty of the fit depends on the 
number of stars included in the fit, and so on the lo- 
cal stellar density, the sky area and range of D used. At 
high latitudes, binning together ten SDSS fields, for stars 
with 10 < D < 19 mag, typical statistical uncertainties 
are 20, 15, 6 and 6 mmag for u — g, g — r, r — z, and 
i — z, respectively. Finding the standard deviations of 
measurements taken near the north Galactic pole, where 
the expected reddening is smooth and small, we find un- 
certainties of 24, 17, 10, and 12 mmag, similar to the fit 
results. 

The presence of dust reddens stars, pushing the blue 
edge redward (Fig. H]). In low-extinction regions where 
all of the dust is closer than about 1 kpc, the observed 
blue tip is uniformly shifted. It is this shift that we 
track in our measurements. In high-extinction regions or 
regions with dust beyond about 1 kpc, we can in principle 
learn about the change in dust column with distance by 
varying D, though we have not done so in this work (see 

The functional form P(x,xq) chosen to fit the blue 
edge of the distribution is usually a close fit only near 
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Figure 3. Example fits for the blue tip of the stellar locus in u — g, g — r, r — i, and i — z for a one degree region around / = 150° and 
b = 20° . A grayscale of the color-magnitude diagram is plotted, with features of the fit superimposed. The diagonal solid lines in the g — r 
plot give the magnitude range within which stars are fitted, corresponding to 10 < D < 19. The histograms give the distribution of stellar 
colors satisfying this cut. The solid error function gives the best fit to this histogram, within the gray dashed vertical lines. The black 
dotted vertical line gives the derived location of the blue tip. The arrow in the upper left of the diagram gives the magnitude and direction 
of the reddening vector in this field, according to SFD. Finally, the right hand axis gives the value of the pdf P; its integral between the 
outer dashed lines is unity. 
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Figure 4. Example fits for the blue tip of the stellar locus in regions of increasing extinction, at I = 150° and b = 80°, 40°, 20° and 10°, 
from left to right. The SFD extinction prediction is indicated by the length of the arrow in the upper left of the diagram. As we approach 
the plane, the stellar density increases. At b = 10°, the brightest stars are not behind all of the dust, and so the color of the blue tip 
becomes redder at fainter magnitudes. 



the blue edge itself, and is rather bad over most of the 
color range. However, experiments with other functional 
forms and attempts to adaptively shrink the range over 
which the fit is performed to a narrow region around the 
edge tend to render the fit less robust. Moreover, be- 
cause P is approximately constant away from the edge, 
the fit results are robust to outliers far from the edge. 
Alternative functional forms systematically shift the lo- 
cation of the edge around slightly (mmags), but not in a 
way that correlates with the dust. 

The blue tip fit can occasionally fail to find the blue 
tip of the stellar locus when too few stars are used. In 
such cases it might identify a single blue star or quasar 
as the blue tip, or latch on to the M-dwarf peak in color 
magnitude diagrams. In the maps we present, with 10 < 
D < 19 and using 10 SDSS fields worth of stars for each 
blue tip fit, fewer than 1% of measurements are affected. 

3.2. Changing Extinction with Distance 

In this work we treat all MSTO stars observed by the 
SDSS with 10 < D < 19 as behind all of the dust. In 
principle, however, by comparing the color of the blue tip 
for nearby stars and distant stars, distant clouds of dust 
could be detected. Such clouds certainly exist, as dust 
has been detected in the halos of other g alaxies and in 
H i clouds outside the dis k in our Galaxy (Me nard et all 
l2010t iWakker et fl/Jfl996h . However, comparing blue tip 
maps made for nearby ranges of D and distant ranges 
of D reveal no readily identifiable structures larger than 
the noise in the halo. We note that if such structures 
vary slowly spatially, then they will be difficult to distin- 
guish from variation in blue tip colors due to changing 



stellar populations. We defer to later work the attempt 
to identify such clouds. 

On most sight lines the great majority of the dust col- 
umn comes from the Galactic disk and is associated with 
the Hi disk, which has a scale height of about 150 pc 
(Kalberla and Kerp 2009). Accordingly, as the brightest 
unsaturated MSTO stars are approximately a kiloparsec 
away, at high Galactic latitudes even the nearest MSTO 
stars observed by the SDSS are behind this dust. 

At low Galactic latitudes, extinction increases rapidly 
and the color of the blue tip of the stellar locus becomes 
dramatically redder as fainter magnitudes are probed 
(Fig. [5j). In these dusty regions, reddened blue tip stars 
may be redder than intrinsically red foreground stars. 
Moreover, blue tip stars in these regions are spread out 
along the reddening vector, making it hard to isolate 
stars of a given reddening and distance using D. We 
attempt to identify and exclude such regions when ana- 
lyzing blue tip maps. 

3.3. Blue Tip Maps 

We measure the blue tip color over the entire SDSS 
footprint, in each color and for a variety of ranges of D 
(Fig. [6]). Here we present maps using blue tip fits to stars 
in a broad range of magnitudes: 10 < D < 19 mag, using 
fits to the stars in 10 adjacent SDSS fields. 

The maps show unmistakable signatures of dust (dark 
clouds in left panels of Fig. [6]) that are almost entirely 
removed when stellar colors are corrected according to 
the SFD map. The residuals at high Galactic latitude 
after extinction correction are dominated by problems 
with the survey calibration — striping in the SDSS scan 
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Figure 5. Example fits for the blue tip of the stellar locus for different magnitude ranges D at / = 150° and b = 10°. The blue tip 
shifts redward as the range of D probed goes fainter, until we run out of stars. Additionally, intrinsically red foreground stars increasingly 
contaminate the blue tip as we go fainter in dusty regions. 
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Figure 6. Blue tip colors over the SDSS footprint, for each SDSS color. The left hand panels give the observed blue tip color, while the 
right hand panels give blue tip colors after correction for dust extinction according to SFD and the coefficients (A a — A^)/(SFD E(B — V)) 
presented in this paper. The upper row of panels gives SFD and the locations of important points in the sky. 

direction and a few runs with bad zero points in u — g blue tip color becomes redder near the Galactic plane, 
(Fig. [7]). There is also a slowly varying residual in which until deep in the Galactic plane when the color shifts 
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very blue. 

The former effect is caused by the various stellar popu- 
lation sampled in different parts of the sky. In the plane, 
the higher metallicity, redder disk population becomes 
increasingly dominant over the halo population. The 
latter effect is a result of the violation of the assump- 
tion that all of the stars are behind all of the dust. Stars 
in front of the dust are dereddened with the full SFD 
dust column, and are therefore rendered extremely blue. 
These stars are then seen by the blue tip algorithm as 
the signature of the blue tip, resulting in a spuriously 
blue color. Even in less dusty regions where at least 
some MSTO stars are behind all of the dust, the method 
can fail if an insufficiently faint range of D is used. In 
this case, MSTO stars behind different amounts of dust 
are grouped together in finding the best fit, blurring the 
blue edge of the stellar locus. Because <jp for the fit is 
fixed, this will result in a bluer than average blue tip 
color, while SFD would track the color of the MSTO 
stars behind all of the dust — the reddest stars in the 
group. These effects, however, are only noticeable in the 
dustiest region of the sky, where \b\ < 15°. 

The Monoceros stream is clearly visible as a feature at 
(a, 5) = (110°, -10°) - (110°, 60°) in the blue tip maps, 
because its MSTO is bluer tha n that of the thick and 
thin disks ([Newberg et a/.ll2QQ2[ ), except in u — g. 

There are two types of survey-systematic induced 
residuals in the maps, both of which are seen as striping 
along the SDSS scan direction. The less important are 
the few, large, ~ 3° wide stripes that are most noticeable 
on the eastern edge of the u — g blue tip maps. These re- 
gions correspond to the footprints of SDSS runs that have 
slight zero point offsets from the rest of the survey. The 
more troubling residual is the universal small-scale strip- 
ing along the scan direction. This corresponds to the 13 
arcminute scale of the SDSS camera columns. We have 
not fully examined this residual, but attribute it to differ- 
ent filter respons es between the 6 SDSS camera columns 
(jDoi et aZ.I [20Toh . The SDSS ubercalibration algorithm 
ties stars of some mean color together, but does not ac- 
count for color terms betwe en the SDSS camera columns 
(jPadmanabhan et a/.ll20M ). As the MSTO stars are sub- 
stantially more blue than the typical SDSS star, they are 
particularly vulnerable to color terms. Complicating this 
explanation of the camera column striping is the obser- 
vation that the best fit camera column offsets derived in 
§4.21 are different for the brightest range of D used than 
for fainter ranges of D, which are mutually compatible. 
We interpret this as a saturation or nonlinearity effect 
that depends on camera column, though we remove any 
stars flagged as saturated from the analysis. These sys- 
tematics need to be addressed in fits to maps of the blue 
tip. 

4. FITTING THE BLUE TIP MAP 

The sensitivity of the blue tip colors to dust redden- 
ing makes for a natural test of the SFD map. We carry 
out the simplest tests of SFD permitted by the blue tip 
maps here, checking the dust map normalization and the 
Ry =3.1 reddening law assumed by SFD for extrapolat- 
ing reddening to colors o ther than B — V. The reddening 
law used i n SFD is fromjO 'Donnelll (jl994f ) . which is sim- 
ilar to the iCardelli et all (j 19891 . CCM) reddening law. 

The SFD dust map is often used as a map of redden- 



ing E(B — V), but is based fundamentally on thermal 
emission from dust. In SFD, the thermal emission is 
converted to reddening based on a single dust map nor- 
malization constant. This constant was derived by com- 
paring the observed reddenings E(B — V) of a sample 
of elliptical galaxies to the SFD te mperature-correcte d 
emission from dust at their locations (jFaber et a/.l ll989). 

We repeat this test, using the colors of the blue tip 
of the stellar locus in place of the B — V color of galax- 
ies. Because of the five band photometry provided by the 
SDSS, we are able to extend the SFD test by addition- 
ally checking the assumed O'Donnell Ry = 3.1 reddening 
law. Because of the number and accuracy of the blue tip 
color measurements, we are able to further look for spa- 
tial variation in the SFD normalization constant and the 
reddening law. The method is ultimately the same as in 
SFD, however: we find the best fit normalization con- 
stants that convert from dust emission (or E(B — V)sfd, 
as these are proportional) to reddening. We extend this 
test to multiple colors and check for variation over the 
sky. 

To perform this test each blue tip map would be fit, 
ideally, as 

(m a - m 6 )biueti P = R a -bE(B - V) S fd + C, (3) 

over a part of the sky, where (m a — ra^biuetip is the blue 
tip map for the color a — 6, R a -b is the normalization 
constant to be measured, and C is the intrinsic blue tip 
color. However, because the blue tip of the stellar locus 
is not a universal color standard, but rather varies with 
position in a way likely substantially covariant with the 
dust map, this fit is impractic al. M oreover, the survey 
striping artifacts discussed in ^3.31 could further throw 
off the fit. To account for these two effects, we instead 
fit 

(m a - m 6 ) b iuetip = Ra-bE(B - V) SF d + Q r (f) + Q (4) 

where Q r (f) is a quadratic for run r in SDSS field num- 
ber / and Ci is the camera column offset for camera 
column i. We fix C\ to be zero to remove the degen- 
eracy with the constant term in Q r (f). The quadratic 
Q r (f) simultaneously accounts for zero point errors in 
the SDSS calibration and slow intrinsic variation in blue 
tip color. 

Because we bin 10 SDSS fields together for each blue 
tip measurement, variations in the dust map on scales 
much smaller than 1.5° in the scan direction will not be 
captured in the blue tip measurements. The blue tip 
method is not a linear operator on the underlying stellar 
colors, so the measured color of the blue tip is not linearly 
related to the mean of E(B — V)sfb in the fields. In this 
work, for each set of fields contributing to a blue tip 
measurement, we use the median E(B — V)sfb in those 
fields as proportional to the expected reddening of that 
measurement. This is not strictly correct, but insofar as 
the dust does not vary too quickly, this filtering seems to 
approximate the correct behavior. 

After finding the best fit R a -b in each color, we re- 
peat the blue tip measurements on the SDSS stars, this 
time with individual stellar colors corrected according to 
the derived fit parameters R a -b, at the full resolution 
of the dust map. This process converges in a few itera- 
tions. The iteration renders the final fits insensitive to 
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Figure 7. Stretched version of Fig. [6] highlighting residuals. The Monoceros stream is prominent at the right hand side of the NGP 
region of the footprint. White to black spans 0.225, 0.11, 0.05, and 0.06 mag. 



the details of the dust map filtering. We have verified the 
insensitivity to the filtering by noting that fits to the un- 
binned blue tip maps using an unfiltered SFD map give 
fit coefficients that are compatible with the fits to the 
binned blue tip maps. 

We model the striping in camera columns by constant 
terms in the fit for each camera column. Empirically 
this correction is satisfactory, but we note that if our 
interpretation of the camera column offsets as derived 
from color terms is correct, the response of each camera 
column to dust is slightly biased. However, we detect 
camera column offsets of ~ 5 mmag over differences in 
stellar colors of a few tenths of magnitudes, implying that 
this bias is less than 1% (Fig. [8]). Furthemore, the bias 
will be different in each camera column and it must be, on 
average, zero. We have further confirmed that our results 
are unaffected by the color terms by separately fitting 
the dust map to each SDSS camera column individually, 
obtaining fits consistent with the global 6 camera column 

fit ( gap . 
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4.1. Fits to Individual SDSS Runs 

In some cases, a single run is long enough and dusty 
enough to individually constrain the dust reddening co- 
efficients R a -b m each color at the 10% level, even when 
fitting for camera column offsets and a quadratic. The 
fits show substantial variation in R a -b (Fig. [9]). Most 
of this variation takes the form of an overall difference 
in dust absorption normalization, rather than as varia- 
tion in the shape of the dust reddening spectrum. Some 
low, outlying values of R u - g suggest that the run is so 
extinguished that substantial numbers of stars are not 
detected in the u band, ruining the fit. The fits remove 
nearly any trace of dust in these runs (Fig. HU|) . 

Fits to the runs individually well constrained by the 
dust map indicate significant variation in normalization 
(~ 20%) as well as, to a lesser extent, variation in dust 
extinction spectrum. Accordingly, in finding a global fit 
to the dust properties, we remove from the fit regions 
that we deem discrepant from the "average" dust prop- 
erties. To do this, first we individually fit each SDSS run. 
Runs with best fit R a -b more than 5cr from the inverse 
variance weighted average of R a -b for all the SDSS runs 
are excluded from the global fit. Additionally, the global 



Figure 8. The camera column offsets used to destripe the SDSS 
blue tip maps, derived from the global blue tip map fits ( ^4.2j) . The 
left-hand axis labels give the range of D to which the offsets apply. 
The right hand labels give the value of the grayscale in mmag. 
The offsets are compatible for different ranges of D, except at the 
brightest magnitudes (10 < D < 15), where nonlinearity effects 
may be coming into play. 

fit is iterated and a-clipped field-by-field at 5<r. 

4.2. Global Fit to All SDSS Runs 

We find the best fit (least squares sense) combination 
of the SFD dust map, constant offsets for each camera 
column other than the first, and quadratics in field num- 
ber for each run to the blue tip map via singular value 
decomposition. We do not fit for the offset for the first 
camera column because changing the offsets for all of the 
camera columns is degenerate with changing the constant 
offsets in the quadratics for each run. Specifically, we find 
the parameters x satisfying 



A T C _1 Ax = A T C _1 b 



(5) 
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Figure 9. R a -b f° r the four SDSS colors, for a selection of indi- 
vidually well constrained SDSS runs. The reddening spectrum has 
a similar shape among the different runs, but there is substantial 
scatter in the overall normalization (left). After forcing the mean 
of Rg-r, R r —ii and Ri—z to match for each run, there is close 
agreement among the various dust extinction spectra (right). Er- 
ror bars account only for the formal statistical uncertainties. The 
thick gray line in each plot gives the global best fit R a -b derived 
in this work. 

which provides the least squares solution to Ax = b. 
Here C is a diagonal covariance matrix, with the vari- 
ances determined by the formal statistical uncertainties 
in the blue tip fit to each field (or binned fields). The 
blue tip color for each field (or binned fields) is in b. The 
design matrix A has dimensions nfi e id x n param . It con- 
tains one column for the SFD dust map values in each 
field, five columns for the offsets for the SDSS camera 
columns two through six, and three columns for each of 
the 792 runs composing the SDSS-III, for the three terms 
in the quadratics for each run. 

The SDSS-III contains imaging data on 1,147,506 
fields. We require that the SDSS scor^3 of each run 
be greater than 0.5 and the PSF FWHM be less than 
1.8 arcseconds in r, which reduces the number of fields 
to 686,554. Excluded runs typically are unphotometric, 
have bad seeing or are Apache Wheel calibration runs. 
The blue tip method is run on all of the remaining fields, 
or binned sets of these fields. It occasionally fails, due 
to an insufficient number of stars or apparent detection 
of the blue tip blueward of —0.75 mag or redward of 2.4, 
2.4, 1.8, and 1.6 mag in u — g, g — r, r — z, and i — z, re- 
spectively. The number of binned fields that successfully 
make it through the blue tip fit depends on the range 
of D used and the color, but in the ranges of D exam- 
ined here it ranges from 684,000 to 677,000. Unbinned 
fields fail for lack of stars at high Galactic latitudes, and 
so fewer fields are successful: in ranges of D of interest, 
between 683,000 and 620,000 fields pass. The median 
number of stars contributing to the fits in each unbinned 
field at 10 < D < 19 is 100. About 10 8 stars enter into 
the fits. 

The global fit produces the coefficients R a -b in Ta- 
ble [1] and at high latitudes leaves few noticeable coher- 
ent residuals (Fig. [TT]) . Nevertheless, a few dust clouds 
clearly are imperfectly subtracted. Most prominently, in 
the north, a few undersubtracted clouds in g — r stand 
out on the eastern side of the north Galactic cap (NGC). 

8 The SDSS score of a field reflects sensitivity to point sources 
and is derived from sky brightness and seeing. It further is capped 
at 0.5 for fields deemed unphotometric or runs fields using binned 
pixels (i.e., Apache Wheel fields). 



Table 1 

R a -b f° r the SDSS colors 



SDSS color 


SFD 


R v = 3.1 O'Donnell 


blue tip 


u-g 


1.362 


1.138 


1.01 ±0.10 


g-r 


1.042 


1.141 


1.01 ±0.08 


r — i 


0.665 


0.616 


0.57 ±0.05 


i — z 


0.607 


0.624 


0.45 ±0.05 



Note. — R a -b according to the original SFD prescrip- 
tion, an Ry =3.1 O'Donnell prescription with updated 
SDSS filter definitions, and the best blue tip global fit co- 
efficients. The formal statistical error bars are negligible; 
those presented here are based on the empirical field-to- 
field variation in best fit extinction law as discussed in J5] 
This gives highly covariant uncertainties because the field- 
to-field variations are dominated by changes in best fit dust 
map normalization (see the covariance matrix in Table [3} • 

Undersubtracted clouds in the northwest and at (a, S) = 
(145°, 0°) cause the runs including them to be excluded 
from the f it. The former f e ature is the largest resid- 
ual in the Peek and Grave maps. Subtraction 
in the south is remarkably clean, though at right ascen- 
sions and declinations of (340°, 20°) and (60°, 0°) there 
are slightly oversubtracted clouds. A cloud at (45°, 20°) 
barely makes it into the SDSS footprint, and is badly 
oversubtracted, causing the exclusion of its run from the 
global fit. In g — r especially, a residual associated with 
the Monoceros stream stands out at the western edge of 
the NGC. We have tried masking this region and repeat- 
ing the fit, but as the Monoceros stream is not correlated 
with the dust, the fit results were unchanged. 

The blue tip color residuals as a function of E(B — 
^)sfd are generally flat (Fig. IT2]). In r-i and i — z, any 
trend with E(B — V) is at less than the 2% level. In u — g 
there is no discernible trend, though the scatter is large. 
In g — r there seem to be disturbing 5-10% trends with 
E(B — V) but they may stabilize around E(B — V) = 
0.4 mag, with reddenings off by only 15 mmag there. 

The blue tip maps have residuals with distributions 
given in Fig. [13l The core of the distribution is Gaus- 
sian, but the wings are non-Gaussian, falling off more 
slowly than would be expected. The blue tip method 
gives median estimated uncertainties as 17.5, 12.5, 5.7, 
and 5.8 mmag, closely reproducing the Gaussian fits to 
the residual distribution, which give 18.1, 12.3, 6.9, and 
7.8 mmag in u — g, g — r, r — z, and i — z, respectively. 

The polynomial terms from the global fit are intended 
to model the slow, intrinsic variation in the blue tip color 
and per-run calibration offsets, and so are of indepen- 
dent interest (Fig. [T4|) . These maps are clearly heavily 
striped. Some of this striping is clearly removing cali- 
bration problems with the SDSS, as in the eastern side 
of the NGC in u — g. The striping at high latitudes in 
the north in other colors is probably substantially an ar- 
tifact of the low signal in these regions, but the typical 
stripe-to-stripe difference there is only ~ 5 mmag, and 
at worst 10 mmag. 

The SDSS imaging contains a number o f runs cross- 
ing t he Galactic plane, as part of SEGUE (jYannv et ali 
2009). All of these runs and a few other low Galactic lat- 
itude runs get identified as discrepant by the fitting al- 
gorithm and discarded from the final fit. In a run- by-run 
treatment, this is inevitable, as at very high extinction 
the blue tip method will fail because the stars are not 



Figure 10. Blue tip colors (left) and residuals (right) for the runs in Fig. [9] in mmag. Coherent signals from the dust are manifest (left). 
After fitting according to the residuals show little structure (right). 



behind all the dust. However, for \b\ > 10°, inspection 
of the blue tip residual maps suggests that the SEGUE 
runs could be well fit also, although we have not tried 
that here. The signal from the dust in the SEGUE runs 
is sufficiently rich to merit treatment separate from the 
bulk of the high Galactic latitude sky, as is done in Berry 
et al (in prep). 

4.3. Fits to different sky regions 

Given the presence of dust-correlated residuals in dif- 
ferent parts of the sky, one wonders whether the fit is 
actually "global" in a useful sense, or is primarily a prod- 
uct of the footprint we have decided to look at. Accord- 
ingly, we have cut the sky into a number of subregions 
and separately found the best fit R a -b in each. To test 
for large scale variations in R a -b, we divide the sky into 
northern (b > 25°) and southern (b < —25°) Galactic 
regions, as well as octants of the sky divided by lines 
of constant Galactic longitude (Fig. [T5]) . To try to test 
variation in dust properties as a function of extinction 
or temperature, we divide the sky into regions of dust of 
different extinctions and temperatures. Finally, to ver- 
ify that we have satisfactorily accounted for the camera 
column offsets, we divide the SDSS survey into its six 
camera columns. 

We restrict to the main global fit region and apply the 
globally derived camera column offsets to the blue tip 
colors. We remove best fit polynomials from the filtered 
SFD dust map and from the blue tip colors, as in 
using the global fit region. We then find the best fit 
of the filtered, polynomial-subtracted dust map to the 
polynomial-subtracted blue tip map, restricted to various 



subsets of the full region. 

The results are encouragingly consistent, given that 
we have already seen that dust reddening normalization 
can vary substantially, along with, to a lesser extent, 
its spectrum (Fig. [9]). We find the dust extinction nor- 
malization and spectrum to be consistent at the 10% 
level. The north/south normalization discrepancy is the 
most surprising, with the north preferring a normaliza- 
tion about 15% larger. We also find that dust with 
0.6 < E(B — V) < 1.0 mag prefers a ~ 15% smaller nor- 
malization than dust with < E(B — V) < 0.05 mag. 
Dust at different temperatures likewise has similar ex- 
tinction normalizations and spectra. 

4.4. Fits to a Mock Star Catalog 

The reliability of the fitting procedure can be tested 
by running the fit on a star catalog from a mock galaxy 
model. The mock catalog was generated with the galfast 
catalog generation code (Juric, in prep). Given a model 
of stellar number density, metallicity, a 3D extinction 
map, the photometric system, and instrumental errors, 
galfast generates realistic mock photometric survey cat- 
alogs. 

As inputs, we used the number density distribution 
parameters from lJuric et all (120081) and the metallic- 
ity distribution from llvezic et all (|2008l ). The three- 
dimensional d ust distribution map was g enerated using 
the model of lAmores and Lepind (|2QQ5h . The gener- 
ated u, g, r, z, and z magnitudes were convolved with 
magnitude-dependent errors representative of SDSS, and 
the final catalog was flux limited at r < 22.5. 

The three-dimensional dust map was modified from 
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Figure 11. Blue tip color residuals after removing best fit linear combination of dust map and polynomials, as described in £14.21 White 
corresponds to +3cr, while black corresponds to — 3cr. 
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Figure 12. Blue tip color residuals with E(B — V)sfd- The 
residuals in u — g, r — i, and i — z reassuringly show no trend 
with E(B — V). In g — r, the residuals have some structure, but 
the positive slope at E(B — V) < 0.2 and the negative one around 
0.2 < E(B — V) < 0.4 are only 5% effects — the residuals are largely 
uncorrelated with E(B — V). 



lAmores and Lepind (J.2005) to contain the small-scale 
clouds seen in SFD. This was performed by scaling 
each line of sight from the model by the ratio E(B — 
V)sfd/ E(B — V) mo dei,iookpc7 so that the model matches 
SFD at 100 kpc. The resulting map contains the ex- 
pected average three dimensional distribution of the dust 
combined with the angular structures present in SFD. 
We construct a mock galaxy using galfast and then 



Figure 13. Fit residuals, in mmag. The core of the distribution 
is well fit by a Gaussian, but the wings fall off more slowly than in a 
Gaussian. The colors r — i and i — z have the tightest distributions, 
likely because the blue tip is less dependent on metallicity in these 
colors. 

construct catalogs of observations of these stars for each 
SDSS run. The resulting catalogs are proccesed by the 
blue tip analysis code in exactly the way that the actual 
SDSS catalogs are processed. As for the actual catalogs, 
the median number of stars per field is ~ 100. 

Blue tip fits performed on the mock galaxy recover the 
R a -b used to within 3%. This verifies that we properly 
account for shifts in the blue tip color due to metallicity. 
However, for these tests we have assumed that SFD cor- 
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Figure 15. R a —b f° r the SDSS colors over different subsets of the SDSS footprint. The first panel gives the coefficients for the north 
(b > 25°) and south (b < —25°). The next panel gives the best fit coefficients for eight slices in Galactic longitude. The third panel gives 
coefficients in regions of sky with increasing extinction. The fourth panel divides the sky into regions of three different temperatures, and 
the fifth panel divides the survey by SDSS camera column. The thick gray line in each panel gives the global best fit R a —b derived in this 
work. 



rectly predicts the line of sight column density and that 
the reddening law is the same everywhere on the sky, at 
least at high Galactic latitudes and in the optical. Inso- 
far as the final fit residuals are largely flat (Fig. [TT]). this 
assumption seems justified. Moreover we can verify that 
the ratios of the R a -b derived independently from SFD 
agree with the ratios of the R a -b derived here. 

4.5. Fitting Ratios of R a -b without SFD 

The preceding analysis has assumed that SFD is a good 
template for the dust. However, we can relax our de- 
pendence on SFD if we restrict our attention to ratios 
between different R a -b- 

In the absence of changing stellar populations and cal- 
ibration errors, blue tip colors will fall along a single line 
given by the reddening vector. By removing a quadratic 
in field number and accounting for camera column off- 
sets as in §4.2| we remove variations in blue tip color not 
associated with reddening. By fitting a line to the result- 
ing polynomial-subtracted blue tip color-color diagrams, 
we measure the ratio of the R a -b- Some striping along 
the SDSS scan direction is evident in the blue tip resid- 
uals despite our attempts to remove it, and this striping 



will be deleterious to the fit results because we cannot 
rely on its being uncorrelated with the dust template we 
fit. Accordingly, we fit lines for each camera column, 
eliminating the effect of striping, and average the results 
(Fig. HB|). 

These color-color fit results are in good agreement with 
the results of £14.21 (Table [2]). The results agree to about 
la with the SFD fit results, where the uncertainties a 
are given by the standard deviations of the fit results 
between the camera columns. The original SFD pre- 
scription is seen to clearly overestimate R u - g /R g - r (see 
£l5.1.3p , while both the original SFD prescription and the 
updated Ry =3.1 O'Donnell reddening law are found to 
underestimate R r -i/Ri- z . 

5. DISCUSSION 

In this section, we use the blue tip fits presented above 
to constrain reddening laws and to test the SFD nor- 
malization and temperature correction. We then briefly 
point out a discrepancy in the color of the dereddened 
blue tip between the north and south, which is plausi- 
bly a consequence of a different stellar population in the 
north than in the south. 




Table 2 

Blue Tip Color-Color Fit Results 



Color ratio 


SFD 


O'Donnell 


blue tip 


color color 


A-u — g/ Ag — r 


1.307 


0.997 


1.005 


1.019 ±0.025 


Ag — rj A r — i 


1.567 


1.852 


1.776 


1.728 ± 0.038 


A r — ij A-i — z 


1.096 


0.987 


1.260 


1.237 ±0.042 



Note. — Blue tip color-color fit results. The results are in 
close agreement with the global blue tip fits, while excluding 
the original SFD reddening law and an Ry =3.1 O'Donnell 
reddening law. These fits do not use SFD as a template and 
so test the shape of the reddening law without reference to 
SFD. 

5.1. Constraining Reddening Laws 

5.1.1. Connecting Reddening Laws to R a -b 

A reddening law gives the extinction A (A) over a range 
of wavelengths A and so makes a prediction for the R a -b 
that we have fit. The extinction in band 6, in the limit 
that the variation in the extinction over that band in 
small, is given by Am^ = A(X e s^,s)E(B — V), with 



A, 



■eff,fc,S 



JdXXS(X)W b (X) 
JdXS(X)W b (X) 



(6) 



Here the system response for band b is given by W&(A), 
S(X) is the source spectrum in photons/s/A, and A(X) 
is the reddening law, normalized to give A(X e s i B,s) — 
^(A e ff,v,s) = 1. The reddening in the color a — b is then 
given by (A(A e ff, ,s) - A(X eSibi s))E(B - V), from which 
it follows that R a - b = (A(A e ff,a,s) - A(X eSibiS ))E(B - 
V)/E(B — V)sfd- Here we have assumed that the vari- 
ation in the extinction over the band pass is small; for 
E(B — V) < 1, this assumption changes the predicted 
R a - b by less than 1% in r — i and i — z, and less than 
5% in g — r and u — g. 

Accordingly, we evaluate the extinctions in each SDSS 
pass band using the SDSS system throughputs and a 
MSTO source spectrum from a Kurucz model with T e ff = 
7QQQK, for a v ariety of reddening laws ( Gu nn et all 1998; 
!Kuruczl[l993h . We examine in particular the CCM, 
O'Donnell, and iFitzpatricl (|1999l F99) reddening laws, 
which are parameterized by Ry = Ay/E(B — V). We 
find the best fit factors Ry to our measured R a - b for 



each of these reddening laws. 

The other free parameter in these fits is the best fit 
normalization N f = E(B - V)/E(B - V)sfd- We do 
not report this parameter directly. The SFD dust map 
is based on a map of thermal emission from dust and 
is correspondingly proportional to Tioo /im, the optical 
depth of dust at 100 /mi. Empirically, the ratio between 
^ioo /xm and E(B — V) is itself a function of Ry, and so 
N' is covariant with Ry. However, the ratio between 
Tl fim and tioo fim depends less on Ry, so we report 

N = A\ jj^m^ predicted/^- 1 /im, SFD- Here A\ predicted 

is the extinction at 1 /im implied by the best fit redden- 
ing law and A\ /xm> sfd is the SFD-predicted extinction, 
extrapolated to 1 /im following the Ry = 3.1 O'Donnell 
reddening law assumed by SFD. The quantity A\ sfd 
is simply 1.32 • E(B - V)sfd. We find that the F99 red- 
dening law provides the best fit to the R a - b we measure, 
so we additionally mention for reference that the F99 red- 
dening law ratio E(B — V)/A\ Mm is 11% larger than the 
O'Donnell ratio E(B — V)/A\ Mm , though this depends 
on the adopted Ry and source spectrum. 

5.1.2. Ry and N for Individual Runs 

We have found that the dust extinction spectrum nor- 
malization and shape vary over the sky (Fig. [9j). We 
can quantify this effect in terms of variation in N and 
Ry by fitting reddening la ws to the runs that individ- 
ually well constrain R a - b ( §4.1)) . We use the F99 red- 
dening law here rather than the traditional CCM or 
O'Donnell reddening law because the latter two redden- 
ing laws (especially the O'Donnell reddening law) predict 
that R r -i ~ Ri-z, which is a poor fit to the data. The 
derived N and Ry vary from run to run, especially as 
a ~ 20% scatter in normalization (Fig. [T7|) . Almost all 
of the runs are consistent with 2.8 < Ry < 3.2. The 
distribution of Ry and N validates our use of 1 /im as 
a refererence wavelength for iV, as the two parameters 
appear uncorrelated. 

5.1.3. R v and N for the Global Fit 

In order to compute best-fit reddening laws for the 
global fit R a -bi we need to account for the covariance 
in the R a - b over the sky in addition to the formal sta- 
tistical uncertainties in the fit. These uncertainties are 
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Figure 17. Best fit Ry and dust map normalization iV for SDSS 
runs that individually well constrain R a —bi examples of which are 
in Fig. [9] Dust map normalizations vary by as much as a factor of 
two, though most are consistent at the 20% level. The large red 
"X" indicates the location of the best fit Ry and normalization 
to the global fit R a -b- Error bars account only for the formal 
statistical uncertainties. 

about a part in one thousand. However, the best fit nor- 
malization can vary at the 10% level from field to field, 
and the part in a thousand uncertainty reflects averaging 
these fluctuations over many fields. While the dust ex- 
tinction spectrum seems relatively constant, it would be 
surprising if the average reddening reflected the redden- 
ing of any particular dust cloud at the part in a thousand 
level, so another method for estimating the uncertainties 
is needed. The uncertainties reported in Table [T] are the 
standard deviations in R a -b fc> r 8 octants in Galactic 
longitude, and do not account for the covariance of the 
measurements in the different colors owing to the chang- 
ing best fit normalization. 

The consequence of adopting uncertainties based on 
the field-to-field variation in the dust extinction spec- 
trum is that the resulting uncertainties on N and Ry 
will describe the range of TV and Ry among the clouds 
within the footprint we analyze. They do not give the 
uncertainties on the best fit N and Ry within our sam- 
ple, which are more tightly constrained. As mentioned, 
the formal statistical uncertainties are about one part in 
a thousand. Uncertainties estimated from using different 
ranges of D for selecting the stars give uncertainties of 
2%, 2%, 1%, and 2% in R a -b for the colors u — g, g — r, 
r — i, i — z, respectively. 

The field- to-field variation in R a -b can be estimated 
from the sample covariance matrix of best fit R a -b in 
different sky regions. We find the best fit R a -b at each 
point in the footprint. Because over most of the foot- 
print the signal-to-noise in the filtered blue tip maps is 
too low to make a reliable determination of R a -b, we in- 
clude nearby pixels in the fit with weights given by a 7° 
smoothed Gaussian. The best fit R a -b is then given by 



Table 3 

Eigenvalues and eigenvectors of the 
reddening covariance matrix 





v u — g 


Vg — r 


V r — i 


Vi-z 


0.155 


0.88 


1.21 


0.45 


0.36 


0.079 


1.23 


-0.51 


-0.41 


-0.79 


0.031 


0.51 


-0.91 


0.61 


1.05 


0.016 


-0.15 


-0.15 


1.35 


-0.83 


Note. 


— Eh 


^envalues 


A and 


eigen- 



vectors v for the sample covariance ma- 
trix. Eigenvectors are normalized so that 
|v| = |R g iobai fit |? to facilitate compari- 
son between R a -b and the first, least well 
constrained eigenvector, which is roughly 
parallel to R a -b- 

where n is the pixel on the sky, i indexes over blue tip 
measurements, D gives the filtered dust map, b and C~ l 
give the blue tip measurement and its inverse variance, 
and Wi^ n is the weight matrix, corresponding to Gaussian 
smoothing with a FWHM of 7°. 

The sample covariance matrix of the resulting maps of 
R a -b is used as the covariance matrix for fits of reddening 
laws to R a -b (Table [3j). Only pixels with estimated un- 
certainty in Ri- Z less than 0.01 are used for computing 
the sample covariance matrix, to avoid overestimating 
the intrinsic variance due to the uncertainty in the esti- 
mates. The largest eigenvalue of the covariance matrix 
(2.4- 10 -2 ) is 94 times larger than the smallest eigenvalue 
(2.6 • 10 -4 ). The eigenvector with the largest eigenvalue 
is roughly parallel to a vector corresponding to chang- 
ing the normalization of the spectrum, and that with the 
smallest eigenvalue corresponds to changing the relative 
amount of reddening in r — i to i — z. This quantitatively 
agrees with our claim that the normalization of the red- 
dening law is poorly constrained relative to its shape. 

With this covariance matrix, the best fit reddening 
laws have Ry = 3.9 ± 0.4, TV = 0.98 ± 0.10 for the 
O'Donnell reddening law, Ry = 3.0±0.4, N = 0.75±0.10 
for the CCM reddening law, and Ry = 3.1 ± 0.2, TV = 
0.78 ± 0.06 for the F99 reddening law (Fig. [18]). The 
X 2 /dof of the O'Donnell, CCM, and F99 reddening laws 
are 7.8, 3.5, and 2.1 respectively. The O'Donnell and 
CCM reddening laws are disfavored because they predict 
smaller reddening differences between r — i and i — z rel- 
ative to the data and the F99 reddening law. The F99 
reddening law fits all of the points reasonably well. 

The O'Donnell reddening law in particular and CCM 
reddening law to a lesser extent provide poor fits to the 
blue tip data, and the best-fit parameters are seriously 
driven by the poorly-matching r — i and i — z data. Ac- 
cordingly, the derived Ry and N for these reddening laws 
are unphysical. 

All of the fits give R u -g substantially less than the 
value given in SFD (Fig. [18|). The values given in the 
SFD appendix were based on preliminary estimates of 
the system response of the SDSS, which varied somewhat 
from the actual sy stem respon s e. Us ing the SDSS sys- 
tem response from lGunn et all ([1998) for the Ry = 3.1 
O'Donnell reddening law gives values different from those 
in SFD by 20% in u — g and 10% in g — r and i — z. 

The predicted R a -b depend somewhat on the source 
spectrum used (Fig. [T9|) . We illustrate the effect of 
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Figure 18. The original SFD reddening coefficients (solid line), 
updated coefficients according to the actual SDSS filters (crosses), 
and the best fit global coefficients found in this work (stars). We 
fit our best fit coefficients with reddening laws according to F99 
(diamonds), CCM (triangles) and O'Donnell (squar es). T he F99 
reddening law seems a good fit to each SDSS band ( ^5.1. 3|) . 
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changing the source spectrum by plotting the derived 
reddening laws for the O'Donn ell, CCM, and F99 red- 
dening laws for Kurucz models ()Kuruczlll993[ ) with T = 
6500, 5500, and 45 00 from the stellar spectra grid of 
iMunari et ali (|2005l ). The blue tip best fit reddening 
law is relatively insensitive to T over this range. 

The maps of R a -b over the sky can be combined to 
make maps of Ry and reddening law normalization over 
the sky (Fig.l20|). To suppress noise in the measurement, 
we have combined the measured values with a prior of 
R v = 3.1 ±0.2 and N = 0.78 ±0.06 to reduce the scatter 
in regions of low signal-to- noise. One feature of the maps 
of Ry and N is that the two maps are substantially un- 
corrected. This again confirms the expectation that the 
ratio of t\ ^m/^ioo does not depend on Ry. 

5.2. Comparison with SFD Normalization 

We can directly compare our predicted E(B — V) with 
E(B — V)sfd by using our best fit reddening law and 
accounting for the difference in source spectrum between 
the galaxies that SFD analyzed and the MSTO stars we 
analyze. We get E(B - V) h iue tip = 0.86 • E(B - F)sfd, 
suggesting SFD overpredicts reddening by about 14%. 
The estimated normalization uncertainty in SFD was 
4%, while we claim 8% uncertainty in our normaliza- 
tion. However, the estimated fractional uncertainty of 
SFD was 16%, and we are now in a position to attribute 
most of this uncertainty to varying best-fit normaliza- 
tion. Accordingly, it is unsurprising that the SFD nor- 
malization differs from ours by 14%, particularly given 
that the footprint we analyze is different from that used 
in SFD, and that we have found north/south normaliza- 
tion differences of about 15%. 

5.3. The SFD Temperature Correction 



Figure 19. O'Donnell, CCM, and F99 reddening laws for dif- 
ferent stars of different temperatures. Ry = 2.6, 3.1, and 3.6 are 
marked with crosses. Ry increases to the left. Ellipses are for 
the best fit blue tip values, with uncertainties computed from the 
covariance matrix of Table [3] 




Figure 20. Maps of Ry (top) and dust map normalization (bot- 
tom) o ver th e SDSS footprint. The map has been combined with a 
prior f 35.1. 3|) to reduce the scatter in regions of low signal-to-noise. 

The final E(B -V) sfd takes the form E(B — V)sfd = 
pI corr X, where p represents a normalization coefficient, 
I CO rr represents the destriped, zodiacal-light subtracted 
IRAS 100 jam flux and X represents a temperature cor- 
rection factor. An error in determining the temperature 
correction factor will lead to dust with a different best 
fit normalization, but with the same reddening spectrum, 
very much like what we see in the blue tip maps. 

Accordingly, the blue tip fit residuals can plausibly be 
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A (mmag) A (mmag) 

Figure 21. Blue tip color residuals versus change in filtered dust 
map E(B — V) from using the FIRAS temperature correction rather 
than the SFD temperature correction (both in mmag). The linear 
trend suggests that the temperature correction may be at fault. 

attributed to errors in the SFD temperature correction 
factor. In §4.3[ we attempted to test the temperature 
correction by verifying that the best fit dust normaliza- 
tion was consistent for regions of different temperature 
according to SFD. However, this procedure assumes that 
the SFD temperature correction adequately distinguishes 
hot and cold dust. At high latitudes, the signal to noise 
of the DIRBE 100 /am and especially the 240 /am map 
are too low to construct a 100 /im/240 fim ratio map 
without substantial filtering. The SFD ratio map was 
constructed by first smoothing the 100 fim and 240 fim 
maps to 1° and then further weighting low S/N pixels 
to a high \b\ average ratio. Insofar as this procedure 
mixes dust of different temperatures together to a single 
reported SFD average temperature, it makes it difficult 
for us to test the SFD temperature correction using the 
SFD temperature map in this way. 

One way to test the accuracy of the temperature cor- 
rection at high \b\ is to compare the SFD tempera- 
ture correction with t he temperature correc tion used for 
the FIRAS dust fits (jFinkbeiner et a/.lll999f ). The SFD 
temperature correction is weighted to constants in the 
north and south in regions of low signal to noise, while 
the FIRAS temperature correction is weighted to a 7° 
smoothed map in these regions. In high signal-to-noise 
regions the two corrections agree. Letting X$fd be the 
SFD temperature correction and Xfiras be the FIRAS 
temperature correction, we can plot blue tip color resid- 
uals versus the change in predicted dust column density 
from switching from the SFD to the FIRAS tempera- 
ture correction: A = E(B - V)sfd(X F t R as/X S fd - 1) 

(Fig. ED- 

Were the FIRAS temperature correction a much better 
predictor of dust temperature than SFD, we would hope 
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Figure 22. The ratio Xfiras/Xsfd in the north (solid) and in 
the south (dashed). The temperature correction is systematically 
higher in the south than in the north by ~ 10%, about compatible 
with the differe nce i n best fit normalization between the north and 
south found in 34.31 

for a linear trend between dust correction and blue tip 
residual, with slope equal to the best fit dust coefficients. 
Instead we find slopes of approximately half the best fit 
coefficients, suggesting that the true temperature map 
is between the SFD and F99 temperature maps. The 
specific values of the coefficients of the linear fits shown 
are not reliable because we have not accounted for the 
uncertainties in A in performing the fit, and so the slopes 
may be underestimated. Regardless, the fact that a clear 
linear trend with significant slope exists provides strong 
evidence that the temperature correction at high \b\ is 
unreliable. 

The FIRAS temperature correction may also explain 
the ~ 15% normalization difference in best fit dust ex- 
tinction s pect rum observed between the north and the 
south in §4.31 The ratio Xfiras/Xsfd in regions with 
E(B - V) < 0.05 is about 10% higher in the south than 
in the north (Fig. [22]). 

5.4. The Dereddened Blue Tip in the North and South 

The color of the blue tip changes because of both dust 
and because of changing stellar populations. If we dered- 
den the blue tip according to SFD and the best fit R a -b of 
this work, we expect the remaining variation in blue tip 
color to be due to changing stellar populations (Fig. [7]). 

A surprising feature of the dereddened blue tip maps 
is that the SGC appears redder than the NGC in g — r, 
r— i, and i — z, especially outside a region in the southeast 
of the SGC (Tabled]). This may be a genuine structure 
in the stellar population in the south, but we are unable 
to distinguish that possibility from calibration errors or 
errors in the dust map. However, because the south tends 
to prefer a smaller SFD normalization than the north, the 
dereddened south would have been expected to be too 
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Table 4 

North/South color asymmetry 



Region 


u- g 


g-r 


r — i 


i — z 


North (mag) 


0.817 


0.228 


0.064 


-0.033 


South (mag) 


0.824 


0.250 


0.071 


-0.020 


Difference (mmag) 


7.6 


21.8 


7.2 


12.4 



Note. — Median blue tip color in the north and 
south, with 40 < |6| < 70. The south is system- 
atically redder than the north. This may be a sign 
of an interesting stellar structure in the south, or a 
calibration problem or SFD error. 



blue rather than too red (Fig. [T5]) . Moreover, the shape 
of the extra reddening in the different SDSS bands does 
not look like any plausible reddening law. 

6. CONCLUSION 

The blue tip of the stellar locus provides a viable color 
standard for testing reddening. It is also a sensitive probe 
of systematics in survey photometry: we detect strip- 
ing artifacts from the SDSS camera columns as well as 
occasional runs with bad zero points. Reddening mea- 
surements based on the stellar locus are limited because 
the stellar locus varies with position on the sky, due to 
changing stellar populations. However, we find that we 
are able to overcome this limitation by looking for small- 
scale fluctuations that correlate with those in the SFD 
dust map. By removing the best fit quadratic from the 
blue tip colors in each SDSS run, both survey systematics 
and sky- varying stellar populations seem adequately ac- 
counted for: Gaussian fits to the blue tip residuals have 
(7 of 18.1, 12.3, 6.9, and 7.8 mmag, compared with the 
median formal statistical uncertainties a of 17.5, 12.5, 
5.7, and 5.8 mmag in u — g, g — r, r — z, and i — z, respec- 
tively. Because our errors are not Gaussian, in terms of 
X 2 /dof we do worse; the fits have x 2 /dof of 1.49, 1.30, 
2.16, and 2.28 in the four colors. 

Using these reddening measurements over the SDSS 
footprint, including, especially, the new, dustier, south- 
ern data, we can sensitively constrain the SFD dust map 
normalization in each SDSS color. The original SFD val- 
ues for Ab/E(B — V) were in error because the filter 
responses used in the SFD appendix did not match the 
eventual filter responses used in the SDSS. After taking 
this into account, our measurements and an Ry = 3.1 
O'Donnell reddening law match closely, except in i — z, 
where the O'Donnell reddening law overpredicts redden- 
ing dramatically. We find an F99 reddening law provides 
a good fit to the data, with Ry = 3.1 and TV = 0.78. We 
recommend the use of this reddening law and normaliza- 
tion for use with SFD and dereddening optical data, or, 
for SDSS data, the use of the constants R a -b presented 
in Table [2 

This result largely vindicates the SFD normalization at 
high Galactic latitudes and low column densities, which 
had been called into question by earlier studies in dustier 
regions. Except in u— g, where the SDSS system response 
assumed by SFD was in error, and in i — z, where the 
O'Donnell reddening law overpredicts reddening, we do 
not find that SFD overpredicts extinction by 40%. In 
g — r, the closest band to the SFD B — V calibration, our 
normalization agrees with the SFD normalization within 
4%. Extrapolating from the SDSS bands to B — V , we 



find that SFD overpredicts E(B — V) by 14%. However, 
because the best fit normalization varies over the sky, 
this result depends on the footprint that we have ana- 
lyzed, and will be different in other areas. The fit results 
indicate that the normalization varies between clouds by 
about 10%. 

We have also made maps of the Ry and N over the 
SDSS footprint. However, we do not yet recommend the 
use of these maps except in regions where the signal- 
to-noise is high, where we use these values to compute 
the R a -b covariance matrix. We are actively investigat- 
ing incorporating similar maps of Ry and N into future 
dust extinction maps, properly combining them with the 
available FIR data. 

The fact that the best fit normalization of the dust map 
varies over the sky points to problems with the SFD tem- 
perature correction. However, the dust extinction spec- 
trum is relatively stable, indicating that at least in re- 
gions with E(B — V) < 0.5, objects can be dereddened in 
the optical assuming a universal extinction law to within 
a few percent accuracy. 

The blue tip colors, dereddened according to the SFD 
dust map with the coefficients from this work, show no 
coherent residuals greater than about 30, 30, 10, and 
10 mmag in the footprint we analyze, which covers most 
of the high-latitude sky. Over much of the sky, residuals 
are within the statistical uncertainties. 

These reddening measurements permit detailed tests of 
dust maps over a large sky area and over a range of dust 
temperatures and column densities. We have been able 
to find clear signs of shortcomings in the SFD tempera- 
ture correction. In future work, we plan to construct new 
dust maps based on dust emission, constrained to best fit 
these measuremen ts. Future surveys like Pan- STARRS 
(Kais er et a/.ll2002l ). the Dark Energy Survey (Flaugher 
| 2QQ5h . and the Large Synoptic Survey Telescope (jTvsonl 
2002) will permit this technique to be used in other bands 
and to fainter magnitudes with better accuracy, allowing 
extended wavelength coverage and admitting dustier re- 
gions near the Galactic plane to be studied. 
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APPENDIX 

USING SFD TO PREDICT EXTINCTION IN SPECIFIC CLOUDS 

Frequently it is useful to estimate the reddening through a particular cloud for which the reddening law is expected 
to be different from the Galactic average, as when the cloud has Ry substantially different from three. The correct 
procedure in such cases is not to use Ay = RyE(B — V)sfd, because in such clouds E(B — V) ^ E(B — V)sfd, as 
SFD is ultimately a map of optical depth at 100 /im, and the ratio of the 100 /im optical depth to E(B — V) is itself 
a function of Ry. 

If the cloud in question has SDSS coverage, the surest footing is to use the measured blue tip colors for that cloud to 
derive the extinction law in that cloud. If SDSS coverage is not available, extrapolation from SFD or some alternative 
method is required. 

In order to extrapolate from SFD, we have found that the best procedure is to take the SFD prediction for E(B — V) 
and transform this to a prediction for extinction at 1 /im. This gives A\ Mm = 0.78/1.32 • E{B — F)sfd, where 0.78 is 
our best fit value for Abiuetip, F99/^-sfd, cd at 1 /im and 1.32 is the ratio E(B — V)/A\ ^ m for an O'Donnell Ry = 3.1 
reddening law. Next, reddening in any other band can be estimated by extrapolating from A\ ^ m according to an F99 
reddening law. Other reddening laws can be used, though they will require constants other than 0.78 in the above 
equation. 
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